Superfluid-Insulator and Roughening Transitions in Domain Walls 
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We have performed quantum Monte Carlo simulations to investigate the superfluid behavior of 
one- and two-dimensional interfaces separating checkerboard solid domains. The system is described 
by the hard-core Bose-Hubbard Hamiltonian with nearest-neighbor interaction. In accordance with 
Ref . [l[ , we find that (i) the interface remains superfluid in a wide range of interaction strength before 
it undergoes a superfluid-insulator transition; (ii) in one dimension, the transition is of the Kosterlitz- 
Thouless type and is accompanied by the roughening transition, driven by proliferation of charge- 
1/2 quasiparticles; (iii) in two dimensions, the transition belongs to the 3D U (1) universality class 
and the interface remains smooth. Similar phenomena are expected for domain walls in quantum 
antiferromagnets . 
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I. INTRODUCTION 



Strongly correlated quantum lattice systems repre- 
sent one of the most exciting and active fields in con- 
densed matter physics. On the fundamental physics side, 
they provide a rich playground to investigate quantum 
phase transitions and study new exotic states of matter 
(for a review see e.g. |3 and as an example of more ex- 
otic systems see e.g. 00]). On the experimental side, 
the possibility of trapping bosons in optical lattices in 
a highly controllable manner makes such systems good 
candidates for applications in a variety of different fields 
such as quantum communication, computing, and preci- 
sion measurements 0, 0, In particular, recently there 
has been a great interest toward studying trapped cold 
polar molecules for which the possibility of controlling 
long-range dipolc-dipole interactions opens up the way 
to realization of novel quantum phases and possible use 
of this system for above mentioned applications (see @ 
and references therein). 

Studies of quantum phases and transitions between 
them are mostly confined to well defined geometries and 
lattice types in a given dimension. It is, however, im- 
portant to recognize that low-dimensional quantum sys- 
tems can also emerge in the form of extended defects in 
a higher dimensional regular structure. Domain walls, 
grain boundaries and dislocations are the most promi- 
nent examples of such systems. Recently, stimulated by 
the observation of non classical moment of inertia in solid 
He-4 9] , superfluid properties of defects in crystals were 
looked at in Refs. [ll, [HI ■ In particular, it has been 
shown that it is possible to have lower dimensional super- 
fluid phases emerging on topologically frustrated defects 
in solid ''He. 

For the scope of this work we are particularly inter- 
ested in the results of Ref. [ij , where the authors stud- 
ied universal properties of the superfluid-insulator (SF-I) 
transition in interfaces separating two checkerboard (CB) 
domains by simulating the classical (d-l-l)-dimensional 
bond-current model [12l | . Here we briefly summarize their 



results. It has been shown that the grain boundary re- 
mains superfluid in a large range of parameters before 
undergoing the SF-I transition. In one dimension, the 
transition is of the Kosterlitz-Thouless (KT) type. An 
argument explaining why one should expect a roughening 
transition to happen simultaneously to the SF-I transi- 
tion has been given: both transitions are driven by prolif- 
eration of one and the same quasiparticle carrying charge 
1/2. In two dimensions, however, the interface remains 
smooth and the transition belongs to the 3D U {!) uni- 
versality class. 

At the qualitative level, the model studied in Ref. [l| 
works for all systems of the same universality class. How- 
ever, it can not be used for making quantitative predic- 
tions regarding realistic quantum bosonic Hamiltonians. 
These predictions are one of the main goals of our pa- 
per. In addition, we present a direct evidence for the fact 
that topological excitations carry charge- 1/2, and quan- 
titatively discuss the connection between roughening and 
SF-I transitions. 

We are interested in studying grain boundaries in 
the CB solid which may be created in the system of hard- 
core cold bosons taken across the SF-CB transition. A 
possible experimental realization is represented by cold 
polar molecules in an optical lattice. In the interesting 
experimental regime, interparticle distances are such that 
double occupancy is strongly suppressed (giving rise to 
hard core bosons) , while the nearest-neighbor interaction 
can be tuned via external electric flelds [1, ■ 

The hard-core extended Bose-Hubbard Hamiltonian 
reads: 
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where a\{ai) is the bosonic creation (annihilation) oper- 
ator, t is the hopping matrix element, V is the nearest 
neighbor repulsion, and /i is the chemical potential. On 
a simple cubic/square lattice the model ([1]) is equiva- 
lent to the spin- 1/2 XXZ antifcrromagnet with exchange 



couplings \ Jx 



\J-u 



2t, — V and magnetic field 
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h — zV/2—fi, where z is the number of nearest neighbors. 
This correspondence makes our work equally relevant for 
studies of domain walls in spin systems, especially in the 
limit of strong anisotropy (or small domain wall width 
when the continuous approximation breaks down). The 
zero magnetic field case, i.e. /i = zV/2, corresponds, in 
bosonic language, to half-integer filling factor. At half 
filling, the ground state of model H]) features the SF and 
CB phases only. The SF phase, characterized by bro- 
ken U{1) symmetry, corresponds to the easy-plane an- 
tiferromagnet with the order parameter = Sx + iSy, 
while the CB order, characterized by broken Z2 symme- 
try, corresponds to the easy-axis antiferromagnet with 
long-range correlations in Sz (here Sx,y,z describe com- 
ponents of the Neel vector). While at a generic filling 
factor the system undergoes phase separation into SF 
and CB phases [11113, at half fiUing the SF-CB transi- 
tion happens at a special higher-symmetry point, V = 2t, 
where the Hamiltonian ([1]) features an 0(3) symmetry. 
In the spin model it corresponds to the Heisenberg point. 
At the Heisenberg point the Neel vector lives on a sphere 
and can point in any direction; it simply rotates from 
equator to the pole when going from the SF to the CB 
phase giving rise to the discontinuous (in terms of order 
parameters) SF-CB transition. 

In what follows we consider a system initially pre- 
pared with two large solid domains with Sz = M and 
Sz = —M respectively, separated by a domain wall. To 
achieve this we consider a system of size Lx = Ly{= 
Lz) = N where N is an even integer, with periodic bound- 
ary conditions (PBC). For the easy-axis Hamiltonian, 
this results in a frustrated CB where two atoms can avoid 
occupying nearest neighbor sites everywhere but in the 
domain wall layer. Our goal is to study properties of 
the resulting (d — l)-dimensional interface embedded in 
a d-dimensional solid using quantum Monte Carlo sim- 
ulations. We focus on the case of zero magnetic field 
(half filling factor), i.e. the chemical potential is fixed to 
/i = zV/2, and show that the interface remains SF (has 
gapless magnons, in magnetic language) well beyond the 
bulk SF-CB transition, for a wide range of V/t. 

Another interesting question concerns the rough- 
ness of the interface (for a brief review of roughen- 
ing transition see e.g. [13 )• To be specific, an inter- 
face (viewed on large scale as a membrane) is said to 
be smooth if its mean square relative displacement, see 
Eq. ^ below, is finite in the thermodynamic limit. On 
the contrary, if the latter diverges with the system size, 
an interface is said to be rough. In this work we present 
numerical results, taken across the SF-I transition, of 
the mean square relative displacement for several system 
sizes. 

The paper is organized as follows. In Subsection II 
A we present results for the ID interface and show that 
the superfluid and roughening transitions happen simul- 
taneously. In Subsection II B we consider a 2D interface 
and calculate the critical point for the SF-I transition. 
We briefly summarize our results in the Conclusion. 



II. RESULTS 

Our simulations are based on the Worm algorithm 
path integral approach [T7||, which allows efficient sam- 
pling of the many-body path winding numbers {W^') in 
imaginary time and space directions. In the path integral 
representation winding numbers have a simple geometri- 
cal meaning: Wt = N is the number of particles in the 
system, and Wx^y,z describe how many times the many- 
body trajectory loops around the system with periodic 
boundary conditions in directions x, y or z. Thus, wind- 
ing number fluctuations in temporal and spatial direc- 
tions define the compressibility and superfluid stiffness, 
respectively. Since we consider the easy-axis Hamiltonian 
a.t V > 2t and half integer filling, i.e. the system is in the 
insulating CB phase, the superfiuidity is restricted to the 
plane of the grain boundary. As a consequence, we expect 
a nonzero value for spatial winding numbers only along 
directions parallel to the grain boundary. In the follow- 
ing we show that the superfluid response persists up to 
a critical value {V/t)c, at which the interface undergoes 
the SF-I transition. In ID interfaces this transition is of 
the Kosterlitz-Thouless type [ll] and is accompanied by 
roughening. In 2D the transition belongs to the 3D U(l) 
universality class and the interface remains smooth. 

In order to calculate critical points (Fig.[2]and Fig. [7] 
below) we have performed simulations in a square system, 
i.e. Lr « Ly^z — L, where L^- = c/S. Here c = Ps/k 
is the sound velocity. In practice, for each system size 
we choose f3 such that (W^) ~ (^Wy z)- This require- 
ment implies c « L/P, with the inverse temperature /3 
scaling with the system size. The last expression for c 
comes from the dependence on winding numbers of the 

■, where d is the system 



superfluid stiffness ps — 
dimensionahty and {W^) — J2i=i d (^i)^ compress- 
ibility K — , where V is the volume. 



One dimensional interface 



Here we consider a system with Lx 



N 



1, 



Ly = iV, which implies a grain boundary along the y- 
direction. Wc start our study at V/t = 2.3, slightly above 
the Heisenberg point V/t = 2. We recall that winding 
number fluctuations in the superfluid phase are described 
by the gaussian distribution (see e.g. which in d = 1 

takes the form P{Wy ) cx eyi-p{—LWy/2(3ps). Similarly 
for winding numbers in imaginary time direction (parti- 
cle number fluctuations) P{W^) cx exp(— /3W^/2Lk). In 
d — 1 the distributions are essentially discrete and the 
best way to extract superfluid stiffness and compressibil- 
ity from the simulated distributions is by considering the 
following ratios: 
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FIG. 1: (Color online) Superfluid stiffness (ps) vs interaction 
strength (V/t) for system sizes Ly = 20 (circles) and Ly — 40 
(triangles) at zero temperature. The dashed line is to guide 
an eye, and the solid line represents ps in the thermodynamic 
limit, from finite size scaling. Error bars are within the symbol 
size. 



FIG. 2; (Color online) Luttinger liquid parameter g as a 
function of V/t. (Error bars are within the symbol size). 
Dashed lines are to guide an eye; the solid line is the re- 
sult of extrapolation to the infinite system size based on the 
Kosterlitz-Thouless theory. Simulation results are consistent 
with the universal jump at = 0.5. The critical value is at 
{V/t)c = 3.57 ± 0.03. 



and 
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We find that — in our system sizes, ensur- 
ing that non-zero values for Wy are due to grain bound- 
ary only. In Fig. [T] we plot the superfluid stiffness as a 
function of the interaction strength V/t, for system sizes 
Ly — 20 (circles) and Ly — 40 (triangles). The temper- 
ature has been chosen much smaller than the finite size 
energy gap, so that the system is effectively at zero tem- 
perature. The grain boundary remains superfluid, and, 
within statistical error bars, we do not see any size effects 
up to V/t > 3.5. Finite-size effects indicate that we are 
approaching the SF-I transition in the interface. 

The ID interface forms a Luttinger liquid and the 
quantum phase transition to the insulating state is of the 
Kosterlitz-Thouless type, characterized by the universal 
jump dX g — gc — where g = ir^psK is the dimen- 

sionless Luttinger hquid parameter (see e.g. [l^l) and 
1 /m is the fractional filling factor. In Fig. [51 we show g as 
a function of V/t for various system sizes (data are taken 
in square systems). Clearly, simulation results are fully 
consistent with the universal jump at gc = 0.5, which im- 
plies an effective filling factor 1/2 in the grain boundary. 
The only logical explanation for fractional filling factor 
when translation symmetry in the bulk is broken with 
doubling of the unit cell volume, is to assume that the 
SF boundary is rough. Rough interface effectively av- 
erages the bulk potential and thus retains the original 
translation symmetry of the lattice. At the microscopic 
level, superfiuidity and roughening are both linked to the 



proliferation of spinon excitations, and we provide addi- 
tional evidence for this explanation below. In the figure, 
the solid line is the result of finite-size scaling based on 
the Kosterlitz-Thouless renormalization group (RG) fiow 
[ri |. The integral form of the equations reads: 



g{L2)/ga 



dt 



= 41n(Li/L2), (4) 



where ^ is a system size independent parameter. Using 
numerical data for g{L) the integral can be solved numer- 
ically and ^{V/t) can be evaluated for each pair of sizes 
(Li,L2)- In the critical region one expects ^ to be a lin- 
ear function of the interaction potential. At the critical 
point, Eq.(|4]) is satisfied by ^ = 1. 

In Fig. [3] we show the solution of Eq. ^ for 
three different pairs of (Li,L2)- The data show a good 
collapse in the vicinity of the critical point {V/t)c — 
(3.57 ± 0.03). The L ^ cx) hmit of Eq. g]) yields 
{g{V/t)/gc){\n[g{V/t)/gc) - i{V/t)) - -1; its solution 
determines the thermodynamic value of the Luttinger liq- 
uid parameter g indicated by the solid line in Figs. [5] and 

m 

We can explicitly verify that topological excitations 
driving the transition in ID carry an effective charge of 
1/2. As it has been argued by Burovski et ai, the propa- 
gation of kinks, or spinous, is achieved by shifting parti- 
cles along the grain boundary. A single particle hopping 
event moves the grain boundary in the transverse direc- 
tion and shifts the kink by two lattice steps. It means 
that in the presence of a gauge field a kink going around 
the system will accumulate half the gauge phase an or- 
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FIG. 3: (Color online) Parameter ^ for different pairs of sys- 
tem sizes {L\,L2) as a function of (V/t). The solid line 
is a linear fit. The transition point is determined from the 
^{V/t) = 1 condition. 



dinary particle will, i.e. its effective quasiparticle charge 
is 1/2. 

In order to show that this is the case we nieasure 
how winding number fluctuations develop in imaginary 
time. The measurement is done for the insulating do- 
main wall when winding number fluctuations are rare 
and one can study individual events. More specifically, 
we monitor Wy and pick configurations with Wy — 1. 
To suppress noise originating from quantum fluctuations 
in the sofid bulk, we "filter" the trajectory by erasing 
hopping events which connect the same nearest neigh- 
bor sites and follow each other in imaginary time (see 
Fig. HI). 

Figure [5] shows an example of the winding num- 
ber trajectory for a system with L — 20 and V/t = 4.2. 
The first step describes system transition to an equivalent 
ground state obtained by shifting interface particles once, 
i.e. with the domain wall shifted by one unit length in 
the transverse direction. At this point the system has two 
choices: going back to the initial state or making another 
transition in the same direction and completing the full 
winding number (which has to be integer). This latter 
case is shown in Fig. [5] because we consider a configura- 
tion with Wy = 1. As one moves closer to the transition 
point spinon excitations become more frequent and start 
overlapping in imaginary time making the whole picture 
less transparent for analysis. 

Next we would like to discuss the connection be- 
tween SF and roughening more quantitatively. In order 
to do so we have calculated the mean square displacement 
(Aa;^) of the grain boundary profile, where the average 
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FIG. 4: Here we show an example on how we "filter" the 
original trajectory a) by erasing hopping events, happening 
on a time scale Ar < 1/V , which connect the same nearest 
neighbor sites. Trajectories b) and c) show the result of 'first 
and second filtering stage' respectively. We stop the process 
at the second stage. 
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FIG. 5: The development of the non-zero winding number 
along the grain boundary in imaginary time. The data refer 
to a system of size Ly = 20 and V/t — 4.2. The presence 
of steps implies that W fluctuations happen in the form of 
instantons, i.e. the transition time from one ground state 
(degenerate in the thermodynamic limit) to another is a rare 
and fast virtual event. The step amplitude ~ 1/2 proves the 
fact that relevant elementary excitations carry charge 1/2. 

is done in both, imaginary time and space: 

{Ax') = j\x^iy,T)dTdy^ . (5) 

Here Ax{y,T) is the instantaneous position of the center 
of the wall (recall that the interface is along the y di- 
rection), calculated from the center of the grain bound- 
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FIG. 6: (Color online) The mean square displacement (mea- 
sured in units of lattice spacing) of the grain boundary profile 
as a function of (V/t). In the insulating phase the data for 
different system sizes collapse, as expected for a smooth inter- 
face. Data start splitting in the vicinity of the critical point. 
In the superfiuid state the interface is rough and (Ax^) de- 
pends on the system size. 



ary at {PLy)^^ J J x{y,T)dTdy. To determine x{y,T), 
we calculate the difference between the "instantaneous" 
densities of two consecutive sites along the x direction, 
dn{i,T) — |n(i,r) — n{i -\- x,t)\. Note that by "instan- 
taneous" density, n{i,T), we mean an average over some 
time interval At such that 1/V <C Ar <C f3, in order 
to eliminate zero point fluctuations. Within the checker- 
board solid, Sn{i, r) ~ 1. In the grain boundary, instead, 
one can have consecutive sites either both occupied or 
empty (with a fractional value oin(i, t) in the supcrfluid 
state of the boundary). Practically, for any given y, we 
start from the first lattice site and determine x{y, t) from 
min((5n(i, r)). 

By definition, (Ax^) is expected to be ~ 1 and sys- 
tem size independent for a smooth interface. For a rough 
interface, however, (Ax^) diverges with the system size. 
Our results are consistent with this expectation as Fig. [6] 
shows. Beyond the transition to the insulating state, re- 
sults for different system sizes overlap within statistical 
errors. On the superfiuid side, instead, (Aa;^) increases 
with the system size L. 



B. Two dimensional interface 

Here we present results referring to a 3D system, i.e. 
2D interface, at half integer filling factor. In this case the 
phase transition belongs to the U{1) universality class in 
three dimensions, indicating a SF — Mott insulator tran- 
sition at integer filling factor. As discussed in Ref. [1], the 
bulk acts as an effective periodic potential which doubles 
the primitive cell in the interface (if it remains smooth 



FIG. 7: (Color online) Finite size scaling for the superfiuid 
stiffness of a 2D interface. From the intersection of curves we 
estimate {V/t)c ~ 2.683 ± 0.003. Solid lines are to guide an 
eye. 



at the transition), bringing the filling factor from 1/2 to 
1. In Fig. [7] we show the finite size scaling of the super- 
fluid stiffness obtained from ps = (Wy + W^) / (3 (recall 
that the transverse direction is x), with /3 scaling with 
the system size. From the intersection of scaled curves 
referring to different system sizes we obtain the critical 
point at {V/t)c = 2.683 ± 0.003. 

In 3D, we did not observe quasiparticles carrying 
fractional charge 1/2. Hopping of single particles along 
the interface remains a local fluctuation; to shift the 
wall one has to create a macroscopic line defect ( "atomic 
step") which is energetically expensive. Since the energy 
barrier between equivalent ground states increases with 
system size the interface remains smooth at T = 0. 



III. CONCLUSIONS 

To summarize, we studied superfluid-insulator tran- 
sitions in 2D and ID domain walls in the bosonic checker- 
board solid. In both cases domain walls remain super- 
fluid well past the bulk SF-CB transition. In ID the 
SF-I transition in the wall is intrinsically linked to the 
interface roughness since both properties are due to the 
proliferation of fractionally charged spinon excitations. 

The minimal description of the system is given by 
the hard-core extended Bose-Hubbard Hamiltonian. Be- 
cause on the bi-partite lattice the latter can be exactly 
mapped onto a spin 1/2 XXZ antiferromagnet model, the 
results presented here are also relevant to spin systems 
(for other lattice types the equivalent spin model is fer- 
romagnetic in the a;j/-plane and antiferromagnetic along 
the z-axis). One can imagine creating domain walls in 
the system of ultra cold bosons in the process of fast 
ramping of the optical potential with several solid seeds 
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nucleated at different locations in the trap. Clearly, the 
dynamics of grain boundaries and sample "annealing" 
will crucially depend on the the SF/roughening transi- 
tions. It might be also possible to observe signatures 
of lower-dimensional coherence in absorption images and 



interference experiments where none are expected for the 
ideal insulating bulk state. 
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